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•  Batteries  in  electric  vehicles  should  be  operated  within  a  narrow  temperature  range. 

•  Tempering  requires  electric  energy  that  has  to  be  provided  by  the  battery. 

•  Optimization  managing  both  demands  is  conducted  using  Pontryagin’s  maximum  principle. 

•  A  causal  tempering  strategy  is  derived  from  the  optimization  results. 

•  This  tempering  strategy  can  be  applied  to  batteries  not  only  in  electric  vehicles. 
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Depending  on  the  actual  battery  temperature,  electrical  power  demands  in  general  have  a  varying  impact 
on  the  life  span  of  a  battery.  As  electrical  energy  provided  by  the  battery  is  needed  to  temper  it,  the 
question  arises  at  which  temperature  which  amount  of  energy  optimally  should  be  utilized  for 
tempering.  Therefore,  the  objective  function  that  has  to  be  optimized  contains  both  the  goal  to  maximize 
life  expectancy  and  to  minimize  the  amount  of  energy  used  for  obtaining  the  first  goal.  In  this  paper, 
Pontryagin’s  maximum  principle  is  used  to  derive  a  causal  control  strategy  from  such  an  objective 
function.  The  derivation  of  the  causal  strategy  includes  the  determination  of  major  factors  that  rule  the 
optimal  solution  calculated  with  the  maximum  principle.  The  optimization  is  calculated  offline  on  a 
desktop  computer  for  all  possible  vehicle  parameters  and  major  factors.  For  the  practical  implementation 
in  the  vehicle,  it  is  sufficient  to  have  the  values  of  the  major  factors  determined  only  roughly  in  advance 
and  the  offline  calculation  results  available.  This  feature  sidesteps  the  drawback  of  several  optimization 
strategies  that  require  the  exact  knowledge  of  the  future  power  demand.  The  resulting  strategy’s 
application  is  not  limited  to  batteries  in  electric  vehicles. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

At  the  moment,  extended  research  is  being  conducted  on  fuel 
cell,  hybrid  electric  and  pure  electric  vehicles  [1—3],  In  this  paper, 
the  focus  is  put  on  electric  vehicles. 

Although  historians  still  argue  who  invented  the  very  first 
electric  vehicle,  electric  vehicles  had  been  invented  as  early  as  in 
the  first  half  of  the  19th  century.  A  brief  history  summary  is  given  in 
[4,5].  Despite  its  early  invention  date,  electric  vehicles  were  su¬ 
perseded  by  internal  combustion  engine  vehicles  at  the  beginning 
of  the  20th  century.  The  reasons  for  this  supersession  were  mainly 
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technical  ones:  the  short  range  and  low  top  speed  of  the  electric 
vehicle  which  were  linked  to  the  poor  battery  storage  capacity.  Due 
to  slow  progress  in  battery  technology  and  high  battery  prices 
resulting  in  high  vehicle  prices,  internal  combustion  engine  vehi¬ 
cles  outperformed  electric  vehicles  [5], 

Nowadays,  the  limited  availability  of  fossil  fuels  and  the 
greenhouse  effect  make  it  necessary  to  further  investigate  low- 
emission  or  preferably  emission-free  propulsion  systems  that  do 
not  depend  on  fossil  fuels.  This  leads  to  a  revival  of  electric  vehicle 
technology  and  extended  research  in  this  field.  Still,  the  factors 
hindering  the  widespread  use  of  electric  vehicles  are  mainly  the 
same  as  back  then  and  are  related  to  the  battery  characteristics  and 
performance:  even  lithium-ion  batteries,  considered  the  best 
available  batteries  on  the  market,  feature  only  low  energy  density, 
leading  to  short  range  and  high  weight  [2, 4, 6, 7],  Additionally,  they 
are  rather  expensive  [8,9],  All  these  characteristics  even  of  the  best 
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available  batteries  make  it  necessary  to  use  an  intelligent  battery 
thermal  and  energy  management  system  [10],  Among  other  tasks 
[10],  battery  management  systems  in  general  contribute  to  opti¬ 
mally  using  the  stored  electric  energy  and  preserving  the  battery, 
increasing  its  life  expectancy.  The  thrifty  use  of  the  stored  energy 
makes  it  possible  to  decrease  the  required  battery  capacity  and/or 
increase  the  vehicle  range,  whereas  increasing  the  life  expectancy 
ensures  long  life  of  the  investment. 

As  electric  vehicles  have  an  electric  machine  and  a  battery, 
regenerative  braking  is  possible.  This  fact  enables  electric  vehicles 
to  come  along  with  less  energy  than,  e.g.,  combustion  engine  ve¬ 
hicles  because  these  cannot  store  and  re-apply  the  energy  gener¬ 
ated  by  the  regenerative  braking  process.  Due  to  the  fact  that  the 
battery  can  be  charged  with  electricity  from  renewable  resources, 
even  CC>2-free  propulsion  is  possible. 

Apart  from  the  electrical  aspect,  the  battery  life  span  also  has  to 
be  considered.  Major  factors  influencing  the  rate  of  battery  aging 
are  temperature,  depth  of  discharge  and  discharge  rate.  [11]  con¬ 
siders  these  factors  for  LiFePCU  cells.  The  discharge  rate  will  not  be 
considered  in  this  paper,  because  it  is  defined  by  the  driver’s  de¬ 
mand,  which  should  be  fulfilled.  Time  effects  on  the  battery  during 
storage  are  also  neglected,  as  these  cannot  be  affected  as  well.  The 
same  electric  power  demand  has  a  varying  impact  on  the  battery 
life  expectancy  depending  on  the  actual  battery  temperature  [9,12], 
For  physical  effects  in  the  battery  causing  degradation,  see  e.g.,  [13]. 
Because  of  the  temperature  effects  impacting  battery  aging,  it 
makes  sense  to  heat  or  cool  the  battery  to  operate  it  in  a  gentle 
temperature  range.  An  optimal  strategy  has  to  be  found  to  balance 
both  the  requirements  energy  saving  and  mild  battery  usage.  The 
need  for  such  an  optimal  strategy  motivates  the  use  of  optimization 
techniques.  Optimization  techniques  have  been  applied  to  various 
tasks  in  the  field  of  automotive  technology:  to  hybrid  electric 
vehicle  power  management  (i.e.,  the  power  split  between  com¬ 
bustion  engine/fuel  cell  and  electric  motor)  [14],  component  sizing 
[14]  and  battery  thermal  performance  optimization  in  electric  ve¬ 
hicles  [15],  to  name  a  few.  Examples  of  optimization  tools  are  fuzzy 
logic  [16],  neural  networks  [17],  heuristic  rules  [18]  and  deter¬ 
ministic  methods.  A  very  powerful  deterministic  tool,  finding  the 
global  optimum  with  respect  to  the  used  value  discretization,  is 
dynamic  programming  (DP)  [19,20],  It  provides  the  optimal  input 
needed  to  reach  optimal  system  behavior.  Though  it  is  reasonable 
to  go  for  the  global  optimum,  the  major  disadvantage  of  DP  is  the 
huge  computational  effort  that  has  to  be  taken  to  find  it  [19],  A 
deterministic  alternative  also  capable  of  finding  the  globally 
optimal  solution  is  Pontryagin’s  maximum  principle  (PMP)  [21], 
PMP  (or  its  dual  formulation,  Pontryagin’s  minimum  principle)  has 
been  applied  successfully  to  various  engineering  problems  in 
automotive  technology,  e.g.,  to  the  power  split  problem  in  hybrid 
electric  vehicles  [21—23]  or  to  finding  an  optimal  cold  startup 
strategy  for  plug-in  hybrid  electric  vehicles  [24],  It  provides 
necessary  conditions  for  finding  the  optimum.  In  some  cases,  these 
can  be  even  sufficient.  The  application  of  PMP  is  comparatively 
convenient  due  to  the  fact  that  it  requires  the  solution  of  differ¬ 
ential  equations  rather  than  calculating  all  possible  solutions  as 
with  the  brute-force  DP  approach.  DP  is  able  to  easily  involve  state 
constraints,  but  these  are  not  required  by  the  present  problem. 

Both  DP  and  PMP  require  the  whole  driving  cycle  to  be  known 
before  starting  the  optimization  process.  Section  6  addresses  this 
drawback.  It  identifies  the  main  parameters  of  the  driving  cycle  that 
have  to  be  known  instead  of  its  detailed  variation  in  time.  With 
these  parameters,  a  close-to-optimal,  causal  strategy  derived  from 
the  PMP  optimization  results  is  applicable. 

The  intention  behind  the  present  work  is  to  abstract  detailed 
models  of,  e.g.,  battery  and  battery  tempering  to  system  level, 
where  the  reasonable  use  of  PMP  is  enabled.  For  a  detailed 


description  of  lithium-ion  batteries,  e.g.,  the  spatial  distribution  of 
temperature  and  current  density,  we  refer  to  [25,26],  Neglecting 
specific  effects  results  in  a  lack  of  accuracy,  but  in  turn,  key  factors 
greatly  impacting  the  optimization  can  be  identified.  Fig.  18  at  the 
very  end  of  this  paper  shows  a  comparison  between  PMP  and  DP. 
For  the  PMP  optimization,  it  is  assumed  that  battery  resistance  and 
open-circuit  voltage  are  constant  and  do  not  depend  on  state-of- 
charge  and  battery  temperature,  whereas  during  the  DP  optimi¬ 
zation,  the  dependence  is  involved.  The  outcome  of  the  PMP  opti¬ 
mization  is  mainly  governed  by  the  initial  temperature,  not  by  the 
values  of  the  aforementioned  two  parameters,  see  also  Section  6. 
Both  PMP  and  DP  optimization  results  hardly  differ,  indicating  that 
the  simplifications  can  be  made  without  a  significant  loss  in 
accuracy. 

Starting  from  the  battery’s  system  model  described  in  Section  2, 
it  will  be  shown  that  no  formal  solution  of  the  PMP  optimization 
(Section  3)  applied  to  the  given  problem  could  be  found.  Still,  the 
problem  can  be  solved  numerically,  see  Section  4.  From  the  results 
of  the  numerical  solution,  presented  in  Section  5,  a  generally 
applicable,  causal  solution  is  derived  in  Section  6. 


2.  Battery  system  model 


To  solve  the  optimization  problem,  the  battery  has  to  be 
modeled  electrically  and  thermally.  The  electric  system  model  of 
the  battery  consists  of  a  voltage  source  with  its  internal  resistance 
connected  in  series,  see  Fig.  1  [27],  The  model  is  also  used  in  similar 
system  level  optimization  scenarios,  e.g.,  [9],  When  it  is  evident 
from  the  context,  the  argument  (t)  will  be  omitted  in  the  following. 
The  open-circuit  voltage  Uoc  depends  on  the  current  state-of- 
charge  SOC,  while  the  internal  resistance  f?bat(SOC,i?)  is  a  function 
of  both  SOC  and  the  battery  temperature  r).  The  terminal  voltage 
l/bat  is  calculated  by 

Ubat=  (V-WW  (1) 

The  equation  for  the  battery  current  fbat  is 

,  Ptrac  +  Ptemp  +  Paux 

'bat - f]— - •  W 


where  Ptrac  stands  for  the  demanded  traction  power,  while  Ptemp  is 
the  electric  power  that  is  needed  to  temper  the  battery.  Paux  is  the 
power  required  by  auxiliary  consumer  loads  other  than  the  battery 
tempering  system.  As  both  Ptrac  and  Paux  are  time-dependent  sys¬ 
tem  parameters,  they  can  be  merged  into  Ptrac.  In  this  paper,  Paux  is 
assumed  to  be  zero.  Eqs.  (1)  and  (2)  combine  to 


-  V^oc  -  4Rbat  (Ptrac  +  Ptemp) 


Only  the  negative  root  is  considered  due  to  the  fact  that  an 
increased  power  demand  requires  a  higher  battery  current.  The 
state-of-charge  SOC  is  defined  by  the  ratio  between  currently 
stored  charge  C(t)  and  nominal  capacity  Cnom. 


HO 


Fig.  1.  Battery  equivalent  circuit  diagram. 
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Using  Eq.  (2),  the  derivative  of  the  si 
expressed  by 


soc(t)  = 


(4) 

n  Eq.  (4)  is 


(6) 


(fbatRbat  +  Q-amb  +  Q-/ 

(»oc--y/^c-4(PtrJC+Ptemp)Bbat) 


Q+Qamb 


r  (v(Fr  +  Fa  +  \n 
'  \(v(Fr  +  Fa  +  irr 


The  PMP  then  provides  necessary  conditions  for  the  calculation 
of  the  optimal  system  input  u(t)  that  minimizes  the  objective 
function  (12).  t0  is  the  beginning  of  the  considered  time  interval, 
while  te  is  its  end  point.  Using  the  definition  of  the  Hamiltonian 


(13) 


y/u$c  4R*x(Pmc-  P,rmp)  Uoc 
—  2Rta,Cnra  -• 

Thermally,  the  battery  is  described  as  lumped  mass  with  heat 
capacity  Cth.bat-  This  leads  to  the  following  equation  for  the  battery 
temperature  r): 


(8) 


In  Eq.  (7),  Qamb  is  the  heat  emission  to  the  environment  while 
the  variable  Q  is  the  heat  power  transferred  to  the  coolant. 

The  electric  power  Pte mp  only  accounts  for  the  electric  power 
necessary  to  produce  the  heat  flow  rate  Q.  at  the  battery. 

The  traction  power  Ptrac  for  plain  road  is  calculated  from  the 
given  velocity  v  by  [28] 

^mv2))/^,  propelling 
2 mv 2))-t),  braking. 

Fr  is  the  (velocity  dependent)  resistance  to  rolling.  Fa  the  (ve¬ 
locity  dependent)  aerodynamic  drag,  m  stands  for  the  vehicle  mass, 
t]  is  the  traction  system  efficiency  factor.  It  depends  on  both  torque 
and  rotational  speed.  For  a  given  driving  cycle  providing  the  vehicle 
speed,  the  required  traction  power  is  calculated  offline. 

3.  Pontryagin’s  maximum  principle 

Pontryagin’s  maximum  principle  (PMP)  is  of  outstanding 
importance  in  the  field  of  optimal  control.  It  provides  necessary 
conditions  for  the  calculation  of  the  optimal  system  input  that 
minimizes  a  given  objective  function.  It  has  been  applied  success¬ 
fully  to  the  power-split  problem  between  internal  combustion 
engine  and  electric  motor  in  hybrid  vehicles,  see,  e.g.,  [29], 

In  the  following,  the  optimal  values  for  input,  state  and  costate 
variables  will  be  indicated  with  the  superscript  o.  For  a  detailed 
description  on  the  dual  minimum  principle,  see  [30],  The  PMP  is 
applied  to  a  nonlinear,  time-variant  state  space  model, 


the  necessary  conditions  are 

*°=!|  =  /(x°(t),u°(t),t),  (14a) 

r  -  -g|o,  04b) 

H(x°,  .  u° .  t)  4  max ,  (14c) 

x°(t0)  =  x0,  (14d) 

x°(te)  =  Xe.  (14e) 

The  last  condition  (14e)  is  only  applicable  when  the  end  value  Xe 
is  given  and  fix.  For  a  free  end  value,  the  transversality  condition 

-  -(s)L  (,5) 

has  to  be  fulfilled.  The  vector  is  the  costate  vector;  each  state  has  a 
corresponding  costate.  Eq.  (14c)  is  called  control  equation  and 
yields  the  optimal  input  u°.  As  not  all  values  for  u°  that  are  calcu¬ 
lated  by  the  PMP  are  feasible  (for  example,  one  or  more  inputs  are 
bounded),  S  in  Eq.  (14c)  is  the  subspace  of  Rp  (p  =  dim(u))  which 
contains  the  valid  values  of  u. 

4.  Optimization  using  Pontryagin’s  maximum  principle 

The  following  arguments  will  be  omitted  for  the  rest  of  this 
paper  when  their  presence  is  clear  from  the  context: 


SOC  (t) 

Pfemp  —  Ftemp^Q(f)^ 

Ptrac  =  Ptrac  (f) 

Uoc(SOC(t)) 

\h(t) 

Fbat(SOC(t),  tf(t)) 

*2  =  Mt) 

m 

H  =  H{x,\j/,u,t) 

QambW- 

4.1.  Pontryagin's  maximum  principle  applied  to  the  system  model 


x(t)  =/(x(t),u(t),t),  (9) 

where  x(t)  is  the  state  vector  and  u(t)  the  input  vector.  It  uses  the 
boundary  conditions 

x(t0)  =  Xq,  (10) 


As  mentioned  before,  the  goal  of  the  optimization  is  to  find  a 
strategy  which  operates  the  battery  preferably  within  a  tempera¬ 
ture  range  in  which  it  is  stressed  least,  while  at  the  same  time  the 
amount  of  electrical  energy  needed  to  temper  the  battery  appro¬ 
priately  is  minimized.  This  leads  to  the  objective  function 


X(te)  =  Xe 

and  the  objective  function 
j  =  h(x(te),te)+ J  fo{x(t)Mt),mt. 


(ii)  j  =  /(-soc  +  vMo^dt. 


(16) 


In  Eq.  (16),  the  integral  of  -SOC  means  that  the  SOC  at  the  end 
time  has  to  be  as  high  as  possible.  /(<5)  is  a  mathematical  quantity 
(12)  that  describes  the  influence  of  temperature  on  the  battery  life;  the 
higher  its  value,  the  more  the  battery  is  stressed,  wi  is  the 
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weighting  coefficient  between  both  effects.  Its  value  results  from 
the  provided  requirements  on  the  battery  life  span  dependent  on 
the  operational  demands,  i.e.,  while  undergoing  a  certain  driving 
cycle  a  fixed  number  of  times,  the  battery  has  to  last  a  predefined 
time  period.  From  the  objective  function  (16),  the  rate  of  heat  flow 
Q(t)  (system  input)  that  minimizes  its  value  has  to  be  determined. 
Substituting  Eq.  (6)  in  Eq.  (16)  results  in 


Uoc-  ^o2c-4(Ptrac  +  Ptemp)Rbat 


Therefore,  the  Hamiltonian  for  this  optimization  problem  is 


„  Uoc-^&-4(Ptrac  +  Ptemp)Rbat 
H  = - 2CnomPbat - WlA^ 

h  (uoc  -  y/U&  -  4(Ptrac  +  Ptemp)  Rbat) 
2CnomPbat 

/  ( Uoc-  VUo2c-4(Pmc+P.«mp)Rb«')  +Q+Qamb\ 
— 1 - 

+  ■ 


(18) 


4.2.  Application  to  a  simplified  system  model 

To  come  to  a  formal  solution,  the  system  model  will  be  simpli¬ 
fied:  Although  state-of-charge  and  battery  temperature  vary 
greatly  during  time,  battery  resistance  Rbat(SOC,il)  and  open-circuit 
voltage  Uoc(SOC)  are  considered  constant,  neglecting  the  influence 
of  SOC  and  il  over  time.  This  simplification  sacrifices  accuracy  in 
order  to  come  a  formal  solution  at  all.  At  the  beginning  of  the 
calculation,  the  parameters  Rbat  and  Uoc  will  be  assigned  the  correct 
value  from  their  respective  lookup  table,  depending  on  the  initial 
values  of  SOC  and  rl.  This  correct  value  will  be  held  constant  during 
calculation.  The  traction  power  Ptrac  is  also  considered  constant.  The 
relation  for  Ptemp  given  in  Eq.  (19)  is  substituted  by  the  general 
Ptemp  (Q)  =  aQ;  the  corresponding  value  for  a,  ah  or  ac,  will  be 
substituted  later  on.  Using  Taylor  series  approximation  for  the 
battery  current  with  different  accuracy  for  SOC  and  i), 


4>at,SOC 


Rbat(W  +  Ptrac)  aQ  + Ptrac 
-'U$c  ■  Uoc 


and 


(24) 


4>at,i? 


(25) 


Considering  the  fact  that  the  battery  can  be  heated  or  cooled,  the 
relation  between  Q  and  Ptemp  is  defined  linearly  by 


there  results  a  simplified  system  model  with  the  equations  for  SOC 


Ptemp  (O.)  j 


“h-Q 

acQ 


for  Q  >  0 . 
for  Q  <  0  ’ 


(19) 


SOCs 


(Rbat(aQ  +  Ptrac)2  +  U2c(aQ  +  Ptrac)) 
C„omU3c 


(26) 


since  for  both  heating  (Q  >  0)  and  cooling  (Q  <  0),  a  positive 
amount  of  electric  energy  has  to  be  provided  by  the  battery.  It  has 
to  be  made  sure  by  higher-level  systems  that  the  total  requested 
power  can  be  delivered  by  the  battery.  The  exact  shape  of  relation 
(19)  mainly  depends  on  the  used  tempering  system.  Due  to  the 
generic  character  of  this  paper,  however,  it  is  beyond  its  scope  to 
investigate  characteristics  of  different  systems,  and  thus,  linear 
dependence  is  assumed. 

Keeping  in  mind  that  the  investigated  system  has  two  states, 
state-of-charge  SOC(t)  and  battery  temperature  tf(t),  from  Eq.  (14b) 
follows  that  there  exist  two  costates,  namely 


JO  8  H  1 

^  “  8SOc|0’ 

(20) 

,0  SHI 

^2  ~  L’ 

(21) 

For  the  costates,  no  initial  conditions  ex 
end  conditions 

fist,  but  from  Eq.  (15),  the 

=  0, 

(22) 

r2  (te)  =  0. 

(23) 

Even  when  the  relations  such  as  Rbat(SOC,i?)  are  assumed  to  be 
linear,  hence  the  derivatives  of  the  Hamiltonian  are  easy  to  calcu¬ 
late,  no  formal  solution  of  the  optimization  problem  could  be 
found.  The  discussion  of  the  corresponding  formulas  would  go 
beyond  the  scope  of  this  paper. 


and  t>,  which  discards  the  factor  a2RbatQ  ,  as  for  the  considered 
parameter  set,  a2RbatQ  -C  U2C  holds: 

Q(2aRbatPtrac  +  u2c)  +  Qambu2c  +  VSa  c 

Wat^c 

The  term  A(i>)  is  defined  by 


(27) 


A(0)  =  A0(tf  -  0O)2,  (28) 

with  Ao  being  a  constant  factor  and  r)o  the  battery  temperature  at 
which  least  stress  is  applied  to  the  battery.  Eq.  (28)  causes  the 
system  to  change  the  battery  temperature  towards  #o;  the  more  the 
actual  battery  temperature  differs  from  ilo,  the  more  effort  is  taken 
to  change  it.  Hence,  the  simplified  system  model  has  the 
Hamiltonian 

Hs  _  v  ^(°Q+M2+uMa<i+M) 

Cnom  UgC 

/Q(2aRbatPtrac  +  U2C)  +  QambU2c  +  Vtrac\ 

+  H  Qh,batU02c  j 

(Rb at  (aQ  +  Ptrac)  "  +  U2C  (aQ  +  Ptrac)  ) 

GwmUfc 

-  Wj  A0(t?  -  t>0)2- 

(29) 

For  the  costate  t/'j,  Eqs.  (20)  and  (22)  hold.  As  the  parameters  Rbat 
and  Uoc  are  considered  constant,  the  dependency  of  the  Hamilto¬ 
nian  Hs  on  SOC  is  neglected  and  hence 
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i ft(ty  =  o  vt.  (30) 

Solving  the  necessary  condition  for  maximizing  the  Hamiltonian 
as  required  by  Eq.  (14c), 


for  the  input  variable  4  yields  the  relation 


(31) 


nO  CnomU0c^2(2aWtrac  +  l&)  2RbatPtrac  +  U2C 
*  2a2CthbatRbat  2aRbat  ' 

Substituting  relation  (32)  in  the  system  equations  and  using  Eq. 
(21)  results  in  a  coupled  system  of  differential  equations  for  the 
state  i}°  and  the  costate  t/'-j : 


•  o  _  0Cnom  (2aRbatPtrac  +  U2.)  2RbatU2cPtrac  +  U(*c 

~ 2a2C2hbatRbatU0c  2aCtbbatRbatLl2c 

2aRbat(RbatPt2rac  +  l&Pfrac  -  Qamb^c) 
2aQh,batRbatUoc 


(33) 


temperature.  This  is  mathematically  correct  as  the  system  model  is 
given  this  way,  but  physically  does  not  describe  reality  correctly. 
Therefore,  this  solution  of  the  optimization  problem  is  inapplicable. 


4.3.  Numerical  solution 


As  a  formal  solution  could  not  be  found,  the  optimization 
problem  will  be  solved  numerically.  The  starting  point  is  the 
Hamiltonian  of  the  original  system  model,  see  Eq.  (18).  Using  Eq. 
(19),  the  Hamiltonian  for  the  numerical  solution  becomes 


""-(«)  - 


HIPtemp=flhQ 

HlptemP=acQ 


for  Q  >  0 
for  Q.  <  0. 


(35) 


For  A(i>),  Eq.  (28)  will  be  used.  As  HnUm(4)  is  a  function 
composed  out  of  continuous  functions,  it  is  continuous  as  well, 
especially  at  Q  =  0,  because  the  left-hand  and  right-hand  limits 
match.  In  contrast  to  the  Hamiltonian  resulting  from  the  original 
system  model,  Rbat(SOC,  ?1)  and  Uoc(SOC)  are  considered  constant. 
Fig.  18  justifies  this  assumption.  Like  in  Subsection  4.2,  Rbat  and  Uoc 
will  be  assigned  the  values  defined  by  the  initial  state-of-charge 
and  the  initial  battery  temperature  and  held  constant  during 
time.  This  again  leads  to  the  conclusion 


'A  2  =  2w1Ao(t?°  —  i)q). 


\ ft(t)  =  0  Vt. 


(36) 


The  one  for  the  SOC, 


Using  Eq.  (21 ),  substituting  H  with  Hnum,  results  in 


a2C2h  hatU2c  -  C2om^2  (2aRbatPtrac  +  U2C 


(34) 


is  not  coupled  with  the  other  two  equations. 

All  three  equations  can  be  solved  formally.  The  above-stated 
optimization  problem  is  the  one  with  the  least  simplifications  for 
which  a  formal  solution  still  could  be  found.  Unfortunately,  it  suf¬ 
fers  from  some  serious  drawbacks:  The  traction  power  is  consid¬ 
ered  constant,  which  means  that  for  involving  a  driving  cycle,  it  has 
to  be  split  up  in  time  periods  in  which  the  required  traction  power 
is  at  least  nearly  constant.  As  the  costate  i/'2  has  to  become  0  at  the 
end  of  each  period  (see  Eq.  (23)),  while  its  initial  value  is  not 
known,  the  solution  stringing  together  all  time  periods  is  not  equal 
to  the  solution  considering  the  driving  cycle  as  a  whole  with 
changing  traction  power  and  the  costate  only  becoming  0  at  the 
end  of  the  total  driving  cycle.  Additionally,  even  within  one  period 
itself,  the  solution  is  not  correct,  as  Eq.  (31 )  does  not  take  into  ac¬ 
count  that  Q.  is  bounded.  This  fact  results  in  a  solution  in  which  for 
example  a  is  set  to  ac  because  of  the  fact  that  the  actual  battery 
temperature  is  above  &o,  but  the  solution  yields  a  very  high  heating 
Q.  This  Q  indeed  raises  the  temperature  even  further  (which  is  not 
in  the  sense  of  the  objective  function  (16)),  but  as  negative  electric 
power  (ac  <  0)  is  needed,  said  in  words,  produced,  it  charges  the 
battery  which  outperforms  the  negative  impact  of  the  rising 


•=  [h  2w1J0(*p  -  tf0)-  (37) 

For  the  calculation  of  the  states  SOC0  and  tf°,  Eqs.  (6)  and  (7)  will 
be  used.  The  input  Q.  is  determined  using  Eq.  (14c), 

max  (38) 

i.e.,  the  optimal  input  Q.  is  the  particular  value  out  of  the  valid  set 
of  inputs  that  maximizes  the  Hamiltonian.  As  a  real  physical  system 
will  be  optimized,  one  has  to  involve  that  only  limited  tempering 
can  be  provided,  Qmin  <  4  <  Qmax  with  4min  <  0  and  Qmax  >  0. 
Practically,  Q.  is  found  by  investigating  each  interval  Qmin  <  Q.  <  0 
and  0  <  Q  <  4max  individually.  This  includes  checking  if  there  is  a 
maximum  of  the  Hamiltonian  within  the  interval;  if  there  is  not, 
check  which  boundary  point  yields  the  largest  value  for  the 
Hamiltonian.  Finally,  both  intervals  have  to  be  compared.  This 
procedure  of  finding  the  optimal  value  corresponds  to  the  one  used 
in  [31], 

5.  Results  of  the  numerical  solution 

For  the  numerical  solution,  the  values  for  Uoc(SOC(t))  and 
Rbat(SOC(t).  i)(t))  are  taken  from  the  lookup  tables  displayed  in 
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Fig.  2.  The  considered  parameter  values  are  those  of  an  electric 
vehicle  lithium-ion  battery;  the  exact  numerical  values  cannot  be 
released  due  to  confidentiality.  Please  note  that  the  horizontal  lines 
coincide  with  the  zero  baseline  of  the  vertical  axes,  i.e., 
Rbat(SOC(t),tf(t))  and  Lfoc(SOC(t)).  This  means  that  there  is  a  lot  of 
variation  in  the  data,  e.g.,  the  resistance  at  low  temperatures  is  a 
multiple  of  the  one  at  high  temperatures. 

Before  the  optimization,  the  exact  parameter  values  are  deter¬ 
mined  depending  on  the  initial  values  of  SOC  and  if  and  kept 
constant  for  the  rest  of  the  calculation. 

Except  for  the  time-dependent  traction  power  Ptrac(t),  all 
parameter  values  are  constant.  For  a  fixed  PtraC,  the  input  Q  de¬ 
pends  only  on  the  costate  fa\  in  the  Hamiltonian  (29),  SOC  does  not 
appear  and  i?  only  appears  in  a  separate  summand,  independent  of 
Q,  and  is  therefore  disregarded  while  maximizing  the  Hamiltonian. 

For  Ptrac  —  5000  W,  Fig.  3  shows  the  optimal  system  input  Q  as  a 
function  of  fa. 

Although  this  graphic  displays  the  function  for  specific  param¬ 
eter  values,  it  will  look  similar  for  different  values  of,  e.g.,  Ptr ac: 
Starting  from  a  negative  fa,  maximum  cooling  power  will  be 
applied.  At  a  certain  point  fa 1, 4  jumps  to  or  steeply  rises  towards 
Q°  =  0.  At  the  next  point,  fa, .  Q°  jumps  to  or  steeply  rises  towards 
Qmax.  From  this  schematic  shape,  some  conclusions  about  the 
general  battery  tempering  strategy  resulting  from  Pontryagin's 
maximum  principle  can  be  drawn.  As  the  costate  fa  has  to  become 
0  at  the  end  of  the  considered  time  period,  the  battery  will  not  be 
tempered  at  the  end  of  the  period.  Depending  on  the  circum¬ 
stances,  battery  tempering  will  be  applied  during  some  time  at  the 
very  beginning.  Section  6  will  discuss  in  detail  the  circumstances 
that,  among  other  things,  determine  the  duration  of  the  tempering 
period.  The  plausibility  of  the  strategy  can  be  illustrated  by  taking  a 
look  at  the  objective  function  (16).  When  a  certain  amount  of 
electric  energy  for  the  purpose  of  battery  tempering  during  the 
driving  cycle  is  provided,  it  is  reasonable  to  spend  this  energy 
rather  early.  Doing  this  results  in  a  lower  value  for  the  penalty 
function  /(??)  which  will  be  integrated  during  the  remaining  time. 
As  the  amount  of  the  used  electric  energy  is  fixed,  the  value  of 

h  =  J  -SOCdt  (39) 


will  be  the  same,  no  matter  at  which  time  the  energy  is  used, 
whereas  the  value  of 


h  =  J  vM(tf)dt 


(40) 


will  be  low  when  the  energy  is  used  early. 

The  actual  calculation  will  be  carried  out  as  follows :  The  coupled 
system  of  differential  equations  (7),  (37)  containing  fa  and  i)  will  be 
solved  simultaneously  using  the  definition  of  the  input  Q.  ,  see 

Q°(V>  2°)/W 


Qr» m  |  1  - ~ 

1  o  ,  o.ooi  V22  0.002 

V-5/K-1 

Fig.  3.  The  optimal  system  input  Q  as  a  function  of  the  optimal  costate  in¬ 


sertion  4.3.  This  calculation  step  provides  the  solutions  for 
fait)  and  Q°(t).  Substituting  Q°(t)  in  Eq.  (6)  results  in  an  uncoupled 
differential  equation  for  SOC°(t),  which  will  be  solved  in  the  second 
step. 

Concerning  the  coupled  system  of  differential  equations,  the 
question  arises  whether  it  should  be  solved  backwards  or  forward 
in  time.  For  t?°,  the  initial  temperature  h°(t  =  0)  is  known,  but  its 
final  value  is  not  known  before  calculation.  The  costate  fa(t)  only 
has  the  boundary  condition  (23),  fa(te)  =  0  for  its  end  point,  while 
its  start  value  is  not  known.  As  Fig.  3  shows,  fa(t)  in  general  has 
quite  low  values,  therefore  its  first  decimals  should  be  calculated. 
Additionally,  it  has  to  turn  exactly  to  0  at  the  end  of  the  time  period. 
In  contrast,  the  actual  initial  battery  temperature  is  determined  by 
measurement;  due  to  the  accuracy  of  this  measurement,  it  is  not 
necessary  to  provide  good  accuracy  while  calculating  i')°(t). 
Considering  these  facts,  the  system  of  differential  equations  will  be 
solved  backwards  in  time,  starting  at  te.  By  solving  it  multiple  times 
with  different  end  temperatures,  the  correct  initial  temperature 
will  be  approached  iteratively.  Once  the  difference  between 
calculated  and  measured  initial  temperature  is  within  0.01  K,  the 
solution  is  accepted. 

Fig.  4  shows  plots  of  the  relevant  variables  for  a  self-defined 
driving  cycle  with  a  duration  of  250  s.  As  can  be  observed,  the 
costate  turns  to  0  at  the  end  of  the  cycle,  giving  the  input  Q  a  shape 
according  to  the  conclusions  drawn  from  Fig.  3.  In  contrast,  Fig.  5 
displays  the  optimization  results  for  the  same  driving  cycle,  but 
using  only  its  first  50  s.  In  this  case,  the  duration  is  so  short  that  it 
does  not  pay  off  to  temper  the  battery. 

The  fact  that  no  battery  tempering  will  be  applied  when  the 
duration  of  the  driving  cycle  is  too  short  also  shows  that  splitting  up 
one  cycle  in  parts  of  short  duration  and  optimizing  each  part 
separately  is  not  optimal.  While  optimizing  the  whole  duration 
leads  to  battery  tempering  in  the  beginning  (optimal  case),  opti¬ 
mizing  each  part  would  lead  to  Q  =  0  for  the  whole  time. 
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Fig.  4.  Relevant  variables  for  a  driving  cycle  with  a  duration  of  250  s:  The  first  plot 
shows  the  required  traction  power  Pafat).  The  following  three  ones  display  the  course 
of  SOC°(t),  of  the  battery  temperature  1 >°(t)  and  of  the  costate  (t),  respectively,  each 
when  the  optimal  input  Q  (t)  (last  plot)  is  applied. 
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Optimizing  each  interval  separately  and  adding  the  objective 
function  values  results  in  a  higher  value  than  the  objective  function 
provides  for  the  unsplit  driving  cycle  optimization. 

6.  Causal  strategy  derived  from  the  numerical  solution 

As  the  costate  has  to  turn  0  at  the  end  of  the  optimization  in¬ 
terval,  the  whole  driving  cycle  (i.e„  the  traction  power  PtracM)  must 
be  known  in  advance  to  be  able  to  apply  PMP.  In  general,  this  is  not 
the  case.  For  this  reason,  a  causal  strategy  is  wanted  that  comes  close 
to  the  optimal  solution,  but  only  needs  a  few  parameters  that  can  be 
at  least  roughly  determined  in  advance.  Keeping  in  mind  that  only 
the  terminal  condition  ^  (te)  =  0  for  the  costate  is  available,  it  is 
necessary  to  determine  the  initial  value  i /'l(to)  to  find  a  causal 
strategy.  Note  that  the  superscript  changed  from  o  to  s,  as  the  causal 
strategy  will  provide  suboptimal  results.  Having  found  the  initial 
value,  it  is  possible  to  solve  the  system  of  differential  equations  (7), 
(37)  forward  in  time,  or,  more  precisely,  in  real-time.  One  method  for 
finding  the  initial  value  will  be  presented  in  the  following. 

For  different  driving  cycles,  the  optimization  will  be  carried  out 
for  known  conditions  (initial  temperature,  initial  state-of-charge, 
Ptrac(O).  These  driving  cycles  include  the  standard  driving  cycles 
named  in  Table  1,  but  some  self-defined  cycles  as  well.  The  standard 
driving  cycles  cover  both  aggressive  and  gentle  driving.  Additional 
information  about  driving  cycles  can  be  found  in  [32],  Driving  sit¬ 
uations  that  are  not  suitable  will  not  be  taken  into  account,  e.g.,  in 
some  cases,  the  high  battery  resistance  at  low  temperatures  makes 
it  impossible  for  the  battery  to  provide  the  required  power. 

Figs.  6-9  show  the  relation  between  the  numerically  determined 
1^2  (t  =  0)  and  both  initial  temperature  and  initial  state-of-charge 
(black  dots).  For  considered  periods  shorter  than  the  original  cy¬ 
cles,  only  the  stated  duration  is  taken  into  account.  At  the  end  of 
this  duration,  the  costate  turns  0  and  the  optimization  is  finished. 


PtracW/W 


Fig.  5.  Relevant  variables  for  the  first  50  s  of  the  cycle  already  used  in  Fig.  4.  Note  that 
in  case  of  such  a  short  duration  it  is  not  optimal  to  temper  the  battery.  As  in  Fig.  4,  the 
first  plot  shows  the  required  traction  power  PMC(t)-  The  following  three  ones  display 
the  course  of  SOC°(t),  of  the  battery  temperature  d°(t)  and  of  the  costate  l/'gft), 
respectively,  each  when  the  optimal  input  Q  (t)  (last  plot)  is  applied. 


Table  1 

Standard  driving  cycles  [32,33], 
Driving  cycle  Abbreviation 


Federal  Test 
Procedure  72 
Federal  Test 
Procedure  75 
Japanese 
10-15  Mode 
New  European 

City  Cycle 
Supplemental 
FTP  Start 
Control  03 
Supplemental 
FTP  United 
States  06 
TOV  SOD  E-Car 
Cycle 


FTP72 

FTP75 

JPN1015 

NEDC 

NYCC 

SFTPSC03 

SFTPUS06 

TSECC 


Distance/m  Duration/s  Average 

speed/km  hr1 

11,997  1369  31.6 

17,787  1874  34.2 

4165  660  22.7 

11,017  1180  33.6 

1903  598  11.5 

5766  596  34.8 

12,894  596  77.9 

60,000  3600  60 


Additionally,  the  plots  include  the  linear  approximation 


f2(t  =  0)  =  Ocyc+bcyc-^t  =  0)  +  Ccyc  •  SOC(t  =  0).  (41) 

The  index  eye  for  the  parameters  a,  b  and  c  only  will  be  used 
when  sampling  points  in  general  are  named;  for  a  particular  driving 
cycle,  it  will  be  replaced  with  the  cycle  name.  This  highlights  the 
fact  that  the  parameter  values  are  different  for  each  cycle.  Later  on, 
the  index  will  be  left  out  when  no  pre-calculated  sampling  points 
are  labeled,  but  the  coefficients  are  used  to  causally  determine 
f2(t  =  0). 

For  all  driving  cycles  up  to  a  cycle  duration  of  about  300  s,  the 
simulated  points  can  be  approximated  quite  well  with  the 
approximation  (41).  It  can  be  seen  from  the  figures  that  (t  —  0) 
hardly  depends  on  SOC(t  =  0).  For  this  reason,  Ccyc  in  Eq.  (41 )  can  be 
set  to  0.  It  still  will  be  used,  though.  The  almost  linear  dependency 
between  \ =  0)  and  the  initial  temperature  can  be  explained 
with  Eq.  (37).  At  a  fixed  cycle  duration,  the  variation  of  j0§(t)  de¬ 
pends  linearly  on  the  temperature  difference  il(t')  -  i)q.  Because  of 
the  rather  high  battery  heat  capacity  and  the  bounded  input  Q_  , 
battery  temperature  only  changes  slowly,  therefore  ^2(t)  is  nearly 
constant  and  ^(t)  is  almost  linear.  As  \p°  ( 0  has  to  turn  0  at  the  end 
time,  the  relation  between  r)(t  =  0)  and  $j(t  =  0)  is  nearly  linear. 
For  increasing  cycle  durations,  the  approximation  becomes  more 
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Fig.  6.  The  graph  shows  the  numerically  determined  initial  values  ^(t  =  0)  (black 
dots)  depending  on  both  initial  temperature  and  initial  state-of-charge  for  the 
JP1015Mode  driving  cycle  with  a  duration  of  50  s.  The  plot  also  includes  the  linear 
approximation  of  the  dependence  according  to  Eq.  (41). 
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Fig.  8.  The  graph  shows  the  numerically  determined  initial  values  1 l%(t  =  0)  (black 
dots)  depending  on  both  initial  temperature  and  initial  state-of-charge  for  the  NEDC 
with  a  duration  of  200  s.  The  plot  also  includes  the  linear  approximation  of  the 
dependence  according  to  Eq.  (41). 


inaccurate,  as  the  change  in  temperature  becomes  noticeable  and 
i ?(£)  in  Eq.  (37)  is  not  constant  any  more.  With  increasing  cycle 
durations,  the  difference  between  simulated  points  and  approxi¬ 
mation  first  becomes  bigger  at  points  close  to  the  optimal  tem¬ 
perature  ??o  =  22  °C,  then  at  points  further  away  from  this 
temperature  (see  also  Fig.  16).  The  impact  of  the  battery  parameters 


Fig.  9.  The  graph  shows  the  numerically  determined  initial  values  ^(t  =  0)  (black 
dots)  depending  on  both  initial  temperature  and  initial  state-of-charge  for  the  TSECC 
with  a  duration  of  250  s.  The  plot  also  includes  the  linear  approximation  of  the 
dependence  according  to  Eq.  (41). 


Fig.  10.  For  each  considered  cycle  and  cycle  duration,  individual  values  of  the 
approximation  coefficients  Qcyc,  bc yc  and  Ccyc  according  to  Eq.  (41)  result  from  the 
simulation  runs.  These  three  plots  show  all  calculated  values  of  the  coefficients 
depending  of  the  considered  cycle  duration  time  te. 
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Fig.  11.  The  same  calculated  coefficients  Ocyc.  bLyc  and  Ccyc  as  in  Fig.  10  are  shown  here, 
but  this  time  depending  on  the  traction  energy  Enac  and  not  on  the  considered 
duration.  Data  points  having  the  same  color  and  mark  belong  to  simulation  runs  of 
different  driving  cycles,  but  covering  the  same  duration.  The  mark  +  belongs  to  a 
duration  of  30  s,  while  cyan  diamonds  indicate  a  duration  of  300  s,  for  instance.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  article.) 
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Fig.  13.  bcyc  jointly  depending  on  te  (see  Fig.  10)  and  £ttac  (see  Fig.  11 )  (black  dots).  The 
second-order  polynomial  approximation  of  b  according  to  Eq.  (42)  is  also  included. 


1,000  2,000  3,000  4,000  5,000  6,000  7,000  8,000 

T^trac/kJ 


resistance  and  open-circuit  voltage  on  $j(t  =  0)  can  be  neglected 
compared  to  the  effect  of  the  initial  temperature. 

For  each  driving  cycle,  the  approximated  dependency  of 
^§(t  =  0)  on  i')(t  =  0)  and  SOC(t  =  0)  can  be  expressed  by  the  pa¬ 
rameters  Ocyc,  bcyc  and  CcyC.  These  differ  by  the  considered  driving 
cycle.  Also,  they  depend  on  the  cycle  duration,  as  already  indicated 
by  the  plots  in  Figs.  4  and  5,  respectively.  In  these  plots,  the  duration 
of  battery  tempering  depends  on  the  duration  of  the  driving  cycle. 

Considering  Fig.  10,  it  becomes  clear  that  the  cycle  duration  is 
not  the  only  parameter  that  determines  the  values  of  acyC,  bcyc  and 
Ccyc.  One  parameter  that  should  be  investigated  is  the  traction  en¬ 
ergy  (i.e.,  the  integral  of  the  traction  power  over  time),  because  the 
battery  heats  itself  while  providing  the  energy.  This  makes  it 
necessary  to  provide  more  cooling  compared  to  an  unstressed 
battery  at  temperatures  higher  than  t?o,  while  at  temperatures 
below  i?o,  less  heating  is  required.  See  Fig.  11  for  the  dependency  of 
the  parameters  on  the  traction  energy.  Figs.  12—14  summarize  the 
dependency  of  a,  b  and  c  on  cycle  duration  and  traction  energy. 
Each  plot  also  includes  the  polynomial  approximation 


x  =  P00,x  +  PlO.x'fe  +  Pot  ,x '  Etrac  +  P20,x  fe  +  Pit  ,x '  te '  Etrac 
+  P02,x'Etrac. 

(42) 

where  x  stands  for  one  of  the  coefficients  a,  b  or  c.  The  second- 
degree  approximation  is  suitable;  first-degree  approximation  is 
possible,  but  less  precise. 

From  Eqs.  (41)  and  (42),  the  formula  for  the  calculation  of 
i p2(t  —  0)  from  the  initial  temperature  i'){t  =  0),  the  initial  state-of- 
charge  SOC(t  =  0),  the  expected  journey  time  te,exp  and  the  ex¬ 
pected  traction  energy  that  will  have  to  be  provided  during  the  trip 
Etrac,exp  is 


=  0)  =  (poo.a  +  PlO,a '  te.exp  +  Pot  ,a '  Etrac.exp  +  P20,a '  fe,exp 
+  Pit, a '  te.exp  'Etrac.exp  +  P02,a '  Etrac.exp)  +  (Poo.h 
+  Pl0,b'te,exp  +Poi,bEtrac,exp  +  P20,ft '  te.exp 
+  Pll,b't'e,exp'Etrac,exp +Po2,6'Etrac,exp)^(t  =  0) 

+  (P00,c  +  PlO.c '  te.exp  +  Pot  ,c '  Etrac.exp 
+  P20,c '  fe,ex p  +  Pll,c '  te.exp  'Etrac.exp 
+  P02.c'E2acexp)  -SOC(t  =  0) 

(43) 

Fig.  15  compares  the  acausal,  optimal  optimization  with  the 
causal,  but  suboptimal  strategy.  Although  the  costate  does  not 
exactly  become  0  when  the  causal  strategy  is  used,  both  curves  are 
very  similar.  Note  that  the  initial  values  $(t  —  0)  —  -5  °C  and 
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Fig.  15.  Comparison  between  acausal,  optimal  (dashed/dotted  gray  line)  and  causal, 
suboptimal  (dotted  black  line)  optimization  result.  Applying  the  acausal  optimization 
results  in  i $j(t  =  te)  turning  to  0,  while  applying  the  causal  strategy  leads  to  a  nonzero 
end  value  of  ij'l- 


SOC(t  —  0)  —  0.45  as  well  as  the  journey  time  275  s  and  traction 
energy  lie  between  the  sampling  points  that  have  been  used  to 
calculate  the  approximation  coefficients. 

Within  the  domain  r ?(t  =  0)  <=  [-30  °C,  60  °C], 

SOC(t=0)<=  [0.1,1],  te.exp  c  [10  s,  300  s],  Etrac.exp  e  [0  kj,  8200  kj],  Eq. 
(43)  is  a  valid  approximation.  For  values  outside  the  mentioned 
domain,  it  becomes  inaccurate  or  even  infeasible.  Fig.  16  shows  the 
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Fig.  16.  Plot  of  the  numerically  determined  initial  values  $(t  =  0),  just  as  in  Figs.  6-9. 
The  driving  cycle  here  is  NEDC  with  a  duration  of  600  s. 
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Fig.  17.  Optimization  of  a  1200  s  segment  of  a  driving  cycle  (FTP75).  The  black  line 
shows  the  optimization  of  the  total  duration,  while  the  dashed  dark  gray  line  shows 
the  separate  optimization  of  periods  of  length  600  s,  updating  the  parameter  values 
after  each  period,  according  to  the  end  values  of  the  finished  one.  The  dashed/dotted 
light  gray  line  displays  the  optimization  of  periods  of  length  300  s,  optimized  analo¬ 
gously  to  the  600  s  periods. 


Fig.  18.  Optimization  of  the  same  driving  cycle  segment  as  in  Fig.  17.  The  dotted  black 
line  shows  the  optimization  of  the  total  duration  with  dynamic  programming,  whereas 
the  gray  line  shows  the  optimization  with  PMP.  The  same  traction  power  PMC(t)  for 
both  optimization  runs  is  plotted  in  black.  The  dynamic  programming  optimization 
involves  updating  the  values  for  RbatfSOC,  d)  and  UocfSOC)  from  the  lookup  table  at 
every  time  step  (1  s).  This  lookup  table  even  takes  different  values  for  RbatfSOC,  it) 
during  charging  and  discharging  into  account.  As  can  be  seen,  even  for  the  long 
duration  of  1200  s,  the  solutions  hardly  differ.  This  permits  leaving  the  parameter 
values  constant  during  each  PMP  optimization  period. 


Table  2 

Final  values  of  the  objective  function  resulting  from  the  optimization  of  a 
1200  s  driving  cycle. 


1200  s  0.180707 

600  s  0.184136 

300  s  0.195548 


initial  values  of  f°r  the  NEDC  at  a  duration  of  600  s.  In  this  case, 
the  linear  relation  between  the  initial  temperature  and  \p°(t  =  0) 
does  not  hold  any  more.  Still,  it  is  possible  to  use  a  second-order 
approximation  similar  to  Eq.  (42)  instead  of  the  linear  one 
described  in  Eq.  (41 ).  It  would  then  be  necessary  to  determine  6 
coefficients  instead  of  aCyc,  bcyc  and  Ccyc- 

Fig.  17  shows  a  long  driving  cycle  of  1200  s  length.  In  addition  to 
the  acausal  optimization  of  the  whole  cycle,  the  total  duration  was 
split  in  periods  of  300  s  and  600  s,  respectively.  Each  period  was 
optimized  acausally,  starting  with  the  respective  end  values  of  the 
preceding  period.  The  parameter  values  of  Rbat(SOC,  i?)  and  Uoc 
were  updated  in  between  the  periods  as  well.  Table  2  shows  the 
objective  function  values  for  the  three  different  approaches.  As 
mentioned  before,  the  optimization  of  the  total  duration  results  in 
the  lowest,  i.e.,  the  optimal,  value.  Optimizing  periods  of  300  s 
provides  the  highest  value.  Fig.  17  points  out  that  it  is  possible  to 
split  the  total  duration  in  periods  with  a  maximal  length  of  300  s 
which  can  be  optimized  causally  using  Eq.  (43),  which  holds  true 
for  up  to  300  s.  When  the  future  driving  conditions  cannot  be 
determined,  it  is  possible  to  append  such  intervals  of  300  s  and 
optimize  each  of  them  individually  with  the  causal  approach. 
Alternatively,  when  the  future  driving  conditions  are  roughly 
known,  one  can  use  higher  approximation  orders  than  the  ones 
used  in  Eq.  (43)  and  deduce  an  approximation  formula  for 
i//|(t  =  0)  that  holds  for  cycle  durations  longer  than  300  s. 

7.  Conclusion 

In  this  paper,  Pontryagin’s  maximum  principle  (PMP)  was  used 
to  find  the  optimal  strategy  to  temper  the  battery  such  that  it  is 
operated  in  a  gentle  temperature  range,  while  least  electric  energy 
is  used  for  tempering.  It  was  pointed  out  that  it  is  not  possible  to 
come  to  a  formal  solution  of  the  optimization  problem.  The  nu¬ 
merical  solution  that  was  calculated  instead  was  used  to  develop  a 
causally  applicable  strategy.  This  causal  strategy  sidesteps  the 
drawback  of  the  PMP  of  only  providing  acausal  solutions.  It  is  only 
necessary  to  determine  some  general  parameters  of  the  planned 
journey  instead  of  knowing  the  exact  curve  of  Ptrac(f)  in  advance. 
These  parameters  have  been  identified. 

In  Section  6,  a  formula  was  derived  for  a  fixed  parameter  set  that 
calculates  the  initial  value  ^(f  =  0)  of  the  causal  costate  i/^(t) 
from  the  parameters  initial  state-of-charge,  initial  battery  tem¬ 
perature,  expected  journey  time  and  expected  traction  energy 
required  during  the  journey.  First,  the  dependence  on  the  initial 
state-of-charge  and  the  initial  battery  temperature  was  modeled 
linearly.  In  the  next  step,  the  coefficients  of  the  linear  approxima¬ 
tion  were  related  to  the  expected  journey  time  and  the  expected 
traction  energy.  For  different  parameter  sets,  it  might  be  necessary 
to  find  a  more  general  formula 

V4(t  =  0)  =  f{SOC(t  =  0),  tf(t  =  0),  Etrac.exp,  fe,exp,  d), 

(44) 

where  /  is  a  suitable  function  and  d  is  a  vector  that  includes  addi¬ 
tional  parameters  that  might  need  to  be  involved  for  different 
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applications.  However,  the  causal  strategy  derived  from  the  nu¬ 
merical  solution  makes  it  possible  to  draw  some  conclusions  that 
are  generally  valid: 

•  By  defining  the  penalty  function  (28)  that  punishes  unfavor¬ 
able  battery  temperature  and  because  of  the  large  heat  capacity 
of  the  battery,  there  is  a  widely  linear  dependency  of  1^2  (t  =  0) 
on  the  initial  temperature.  This  dependency  covers  the 
considerable  dependency  of  the  battery  resistance  on  the 
initial  temperature,  see  Fig.  2. 

•  The  same  holds  for  the  initial  state-of-charge:  Although  Rbat(- 
SOC,  if)  and  Uac(SOC)  depend  on  SOC  in  the  considered  case 
(see  Fig.  2),  (t  =  0)  hardly  depends  on  SOC(t=  0). 

•  The  expected  journey  time  and  the  expected  traction  energy 
obviously  are  significant  parameters  that  must  be  used  to 
approximately  calculate  =  0)  when  the  detailed  curve  of 
Ftrac(f)  is  not  known.  The  expected  journey  time  substantially 
determines  whether  and  how  long  the  battery  needs  to  be 
tempered.  The  traction  power  provided  by  the  battery  warms 
it,  but  it  is  not  important  to  know  exactly  how  much  power  has 
to  be  provided  at  which  point  in  time. 

It  is  possible  not  to  consider  the  traction  energy  (i.e.,  the  in¬ 
tegral  of  the  traction  power  over  time),  but  the  integral  of  the 
absolute  value  of  the  traction  power  as  parameter  that  needs  to 
be  involved.  This  is  due  to  the  fact  that  the  battery  current  heats 
the  battery,  no  matter  if  it  is  provided  by  the  battery  or  if  it  flows 
into  the  battery  to  store  recuperated  energy.  It  has  to  be 
considered,  though,  that  the  recuperated  current  raises  the  state- 
of-charge,  which  also  impacts  the  objective  function,  just  like  the 
rising  battery  temperature.  As  the  dependency  of  ^(t  =  0)  on 
the  traction  energy  is  nonambiguous,  it  is  reasonable  to  consider 
the  traction  energy  rather  than  the  integrated  absolute  traction 
power. 

The  development  of  Eq.  (44),  or  its  special  case  for  the 
considered  battery,  Eq.  (43),  makes  it  possible  to  apply  a  strategy 
for  battery  tempering  which  during  the  trip  provides  a  value  for 
the  objective  function  that  comes  close  to  the  optimal  one.  The 
strategy  can  be  implemented  as  follows.  If  the  drive  is  already 
known,  e.g.,  if  it  has  been  driven  before  or  the  destination  was  fed 
into  the  route  guidance  system,  expected  journey  time  and  trac¬ 
tion  energy  can  be  determined  quite  accurately.  When  neither  the 
driver  nor  the  route  guidance  system  provides  the  required  in¬ 
formation,  it  can  be  determined  roughly  by  stochastical  models. 
From  the  two  parameters  as  well  as  the  measured  initial  battery 
temperature  and  the  initial  state-of-charge,  ^(t  =  0)  is  calcu¬ 
lated.  The  differential  equation  (or  its  discrete  counterpart,  the 
difference  equation)  for  ^|(t)  is  then  solved  in  real  time,  using  the 
continuously  measured  battery  temperature.  From  the  costate  and 
the  actually  measured  traction  power  Ptrac(t),  the  input  Q  (t)  is 
calculated  and  applied. 
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